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ABSTRACT 

A collection of Johnson/Cousins photometry for stars with known [Fe/H] is used to generate 
color-color relations that include the abundance dependence. Literature temperature and bolomet- 
ric cor rection dependences a re attached to the color relations. The JHK colors are transformed 
to the iBessell fc Brettl l)1988ft homogenized system. The main result of this work is the tabulation 
of seven colors and the U-band bolometric correction as a function of T e ff, log g, and [Fe/H] for 
— 1.06 < V — K < 10.2 and an accompanying interpolation program. Improvements to the present 
calibration would involve filling photometry gaps, obtaining more accurate and on-system photometry, 
knowing better log g and [Fe/H] values, improving the statistics for data-impoverished groups of stars 
such as K dwarfs, applying small tweaks in the processing pipeline, and obtaining better empirical 
temperature and bolometric correction relations, especially for supergiants and M stars. A way to 
estimate dust extinction from M dwarf colors is pointed out. 
Subject headings: stars: fundamental parameters — stars: general 



1. INTRODUCTION 

Stellar evolutionary tracks and isochrones calculate 
physical radius, luminosity, and effective temperature. 
In order to compare with observable quantities, almost 
always magnitudes and colors, a transformation is essen- 
tial. The need for such transformations is also felt when 
integrated light models (population synthesis models) are 
constructed for comparisons to colors from galaxies and 
star clusters. Color-temperature transformations are also 
used in spectral abundance analysis, design of observing 
strategies, computation of selection effects, and a host of 
incidental astronomical problems. 

One way to approach this problem is to calculate line- 
blanketed synthetic spectra and integrate under filter 
transmission functions to get fluxes, which are then ze- 
roed b y co mpariso n with Vega or other convenient stan - 
dard frBuser fc Kuruczl Il992t iBell fc Gustafssonl ["19891. 
This is a convenient approach, but vulnerable to errors 
in all of the steps of the process: incorrect atmosphere 
structures, incorrect or incomplete line lists, incorrect 
treatment of convection, turbulence, non-LTE effects, or 
line broadening, incorrect filter transmissions, inaccurate 
spectrophotometry of the comparison star or stars, and 
finally, the photometry of the comparison star itself. Ex- 
amples of suspicious circumstances include the fact that 
the blue edge of the U filter is set by the earth's atmo- 
sphere and will inevitably change with time and place, 
the fact that half th e lines in the s olar spectrum have 
yet to be identified l|Kuruczl 1199 2a). and the fact that 
the absolute flu x calibration for stars is uncertain by 
about 5% [e.g. iBerriman et al.1 ll 19921)1 Recent exam- 




Despite the need for empirical alternatives, the present 
paper appears to be the first attempt at empirically cali- 
brating many colors for all stellar temperatures as a func- 
tion of abundan ce and gravity. The inspiration for this 
work comes from[Grccn ( 1988|). Green describes a global 
colot-T e ff calibration generated for attachment to the 



Revised Yale Isochrones ijGreen et al. 1 119871 that pro- 
vides colors tabulated for a (long) list of temperatures, 
surface gravities, and [Fe/H] values. The strategy used 
by Green was to begin with empirical color-color rela- 
tions for solar-mctallicity stars, and adopt the ridge line 
as the starting place. Then one attaches a color-T e // re- 
lation and adds [Fe/H] and gravity dependence by work- 
ing differentially within synthetic color tables. The ap- 
proach here is to do the same, except that the gravity and 
abundance dependences are fit to the stars themselves. 
Synthetic colors are used at very low weight to guide the 
fits where there are few or no stars in the sample (e.g. 
metal-poor M giants). Only oxygen-rich stars are con- 
sidered here. Colo r-temperature relation s for carbon-rich 
giants are given in lBergeat et al.1 1)20011) . 

This paper is divided into sections on procedure, liter- 
ature comparisons, and a concluding section. Supporting 
material (color-temperature table and interpolation pro- 
gram) is available at http://astro.wsu.edu/models/ 

2. PROCEDURE 

2.1. Stellar Data 

The nucleus of the photometr y cata log is the 
compilation of iMorel fc Magnenat I i)1978ft . which is 
firmly Johnson-system. Several ot her pho t ometr y 
sources were i ncluded. These include |Veeder I 111 9 741. 
IBessell llT99ll). iStetson I [|1981ft Ida Costa fc Armandrofll 
(1990), iWood. Bessell. fc Fox I ill 983ft . [ Wood fc Besselll 
i!983f>. Ivon Braun et al. I lll99Sft7 iCarnevI J1 9831. 
ICohen et alJ d!978l>. lElias et allf (I1982T). lElias et all 
J1985H. IFrogel et all II1978H. IFrogel et al] lll979ft 



iPersson et al. I lU980ft. ICohen et al.l I l9801. Frogel et al J 
119811). Ida Costa et al] dl981l>. ICohen k FrogeTTl 9821 , 
IFrogel et al] tl 983ft IFrogel et all l|1990ft . iLeggett et alJ 
l|2001ft . lDahnet alJ l|2002ft . and the 2MASS point source 
catalog. 

As a first try at assign ing abundance measurements t o 
the list of 4163 stars, the lCavrel de Strobel et alJ feOOl). 
iMcWilliaml l|199flft . and lEdvardsson et alJ IJ199.1) abun- 
dance catalogs were consulted. M giants with good pho- 
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tometry were artificially assigned an [Fe/H] of zero ex- 
cept those that belong to clusters, in which case the clus- 
ter metallicity was adopted. M dwarfs were assigned 
an [Fe/H] based on the ir kinematics, most of which 
came from iVeeder I (|1974l) . Young disk objects were as- 
signed [Fe/H] = -0.1, old disk -0.5, and halo -1.5. 
Cluster stars naturally inherited the cluster metallicity. 
Cluster abund ances came from mostly sec ondhand com- 
pilations llWorthev et al. lfl994t Ivon Braun et al. lll998t 
iFrogel et alJll978|). An abundance of +0.3 was adopted 
for NGC 6791 llWorthev & Jowett 1 l200l . LMC field 
stars were assigned —0.3 and stars in the SMC —0.6. 
Some supergiants and very hot stars were artificially as- 
signed [Fe/H] = when no abundance was available, but 
many had literature abundances. Unfortunately, a com- 
plete citing of the abundance sources cannot be given, 
as notes on some of the (perhaps 5%) abundance assign- 
ments have been lost. A total of 1954 useable stars had 
abundances, although the number is considerably less for 
any given photometric color. Odd holes appear in the fi- 
nal data set. For inst a nce, t he primary source for M 
dwarf colors is IVeeder I l)1974|) from which J-band data 
is missing. J7-band data is hard to find for cool stars. 
Available i?-band data has gaps as well. 

Solar metall icity m ean relations for all spect ral types 
from Johnson (1966) and Bcsscll & Brett (1988) were in- 
cluded in the list, with [Fe/H] = —0.1 assumed for these 
"stars" . 

2.2. Photometric Systems 

Perhaps the reason no one has done this exercise be- 
fore is because photometric systems are such a mess. The 
iMorel fe Magnenat I IJ1978J) catalog is on the "Johnson" 
photometric system. For colors involving RI, the target 
system was "Cousins" and we applied the tranformation 
equations of Bcsscll (1979) and Bcsscll (1983) to trans- 
form the Johnson data exc ept for R—I, f or which we used 
a tracing of Figure 3 from lBesselll (J198 3) rather than the 
formula given in the paper. Additional optical data that 
was already on the Cousins system was left there. 

Infrared data was imported from 5 different systems 
(and this is mild compared to the number of systems that 
have proliferated over the years). As a target system, w e 
chose the homogenized system of lBessell fc Brettl ()1988|) . 
Transformations from Johnson-system, CIT-system, and 
AAO- system were used as provided in Bcsscll & Brett 
(jl988t) . Some 2MASS data, mostly attached to NGC 
679 1 stars in the prese nt stellar catalog, were tranformed 
via | Carpenterl l)2001j) formulae to the iBessell fc Brettl 
(1988) system. 

Corrections for i nterst ellar extinction were done using 
the Cardelli et al. ( 1989) extinction curve. Note that cor- 
rections are applied differently for Johnson RI than for 
Cousins RI since the filters are at substantially different 
wavelengths. Such wavelength differences cause negligi- 
ble correction differences at infrared wavelengths. 

2.3. Color-color Fitting 

The photometrically-homogenized, dereddened stellar 
data were then presented to a series of additional pro- 
cessing steps. A multivariate polynomial fi tting pro- 
gram (a modification of the one used in Worthev et al. 
( 1994) to fit spectral indices as a function of stellar at- 
mosphere parameters) was applied to the data. The de- 
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Fig. 1. — The U — B, V—K color-color diagram for uncullcd stars. 
Stars have different symbol types for metal rich ([Fe/H] > —0.2), 
metal-poor ([Fe/H] < —1.2), and intermediate abundance ranges. 
Calibrations for typical giant and typical dwarf gravities are drawn 
in solid for [Fe/H] = 0, dashed for [Fe/H] = —1, and dotted for 
[Fe/H] = —2. Most carbon stars (asterisks) are not plotted as they 
stretch beyond the plot limits along a line from the plotted ones 
up to (V -K,U -B) ss (6,6). 



pendent variable chosen was V — K because it is mono- 
tonically increasing with temperature and insensitive to 
abundance. V — K is a fabulous temperature indicator 
in stars cooler than the sun, and, because of its mono- 
tonicity, can still serve as a temperature-like variable for 
hotter stars. Terms of up to order (V — K) 6 could be in- 
cluded, and up to quadratic terms of log g, [Fe/H], and 
cross-terms. Chemically peculiar stars such as carbon 
stars were excluded from the fit. The color range was 
divided into 5 overlapping sections and each range was 
independently fit. The specific polynomial terms allowed 
in the fit could be different for each section. This allowed, 
for instance, [Fe/H] sensitivity to be manually phased out 
if desired. The fits were done many times. Outlier data 
points were rejected manually with the aid of a graphical 
interface that allowed the name and parameters of each 
star to be scrutinized before rejection. 

The final polynomials were combined in tabular form, 
using a weighted-mean scheme wherein the middle of 
each V — K section was weighted strongly compared to 
the edges of each section. [Fe/H] and log g were tab- 
ulated in 0.5 dex intervals, —2.5 < [Fe/H] < 0.5, and 
—0.5 < log g < 5.5. The resultant color-color relations 
are illustrated in Figures H El El H El and El All stars, 
even if they were rejected during the fitting process, are 
included in the figures. Carbon stars are included, but 
only for illustrative purposes; fits were not attempted. 

For very hot stars of O and B spectral type an addi- 
tional color-color table that was crafted by hand from 
either summary color-spectral type relations or from 
our own color-color relations was employed to refine the 
color-color rel ations for the hottest stars. The source s 
consulted werelSchmidt-Kalerl l|1982D ; iTokunaga 1 12000) ; 
iVacca et al.l l|1996l) . and our own color-color plots. There 
is basically no abundance leverage for very hot stars, so 
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Fig. 2. — The B—V, V—K color-color diagram for uncullcd stars. 
Stars have different symbol types for metal rich ([Fe/H] > —0.2), 
metal-poor ([Fe/H] < —1.2), and intermediate abundance ranges. 
Calibrations for typical giant and typical dwarf gravities are drawn 
in solid for [Fe/H] = 0, dashed for [Fe/H] = — 1, and dotted for 
[Fe/H] = —2. Carbon stars are shown as asterisks. 




Fig. 4. — The V — I, V — K color-color diagram for uncullcd 
stars. Symbols and line styles are as in Figure |21 
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Fig. 5. — The J — K, V — K color-color diagram for uncullcd 
stars. Symbols and line styles are as in Figure |21 



Fig. 3. — The V — R, V — K color-color diagram for uncullcd 
stars. Symbols and line styles are as in Figure |21 



we assumed a zero metallicity dependence. These tabu- 
lated average values were included in the polynomial fits 
as if they were individual stars. 

In the hot star regime, two uncertain areas came to 
light that deserve mention as regards d warf vs. super- 
giant colors. First, in U — B, Schmidt-Kaler (1982) 
data imply a large and distinct color separation be- 
tween dwarfs and supergiants, but the (few) supergiants 
available in our list did not follow the literature trend. 
Nor were the polynomials flexible enough to track these 
changes, mostly because, for O stars, the difference in 
surface gravity is very minor [4.15 (dwarfs ) vs. 4.09 (su- 
pergiants) according to iVaccaet al. (199q)]. In the end, 
we performed a weighted average between the polynomial 



fits and the tabulated values. There is probably consid- 
erable uncertainty left in the supergiant U — B colors, 
perhaps several tenths of a magnitude. This is one area 
that could be vastly improved with more photometry, 
with the caveats that reddening is often a huge factor 
for these intrinsically bright, usually distant stars and 
the fact that fast rotation introduces an inclination an- 
gle dependence in the colors. Users wishing to avoid this 
entirely may want to feed our interpolation program arti- 
ficially high gravities for stars hotter than about 9000 K. 
The sec ond area of debate was that the tabulated H — K 
colors of iTokunaga I l)2000D for O supergiants were about 
0.09 mag redder than for O dwarfs. In this case, we saw 
no trace of such a trend in our data: a few stars were that 
red, but they were all dwarfs. We allowed the polynomial 
fit (which, in that regime, was a function of temperature 
alone) to determine the final color-color relation. In the 
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Fig. 6. — The H — K , V — K color-color diagram for uncullcd 
stars. Symbols and line styles are as in Figure 1^1 



middle of the temperature range, a small gravity depen- 
dence was indicated, but no dependence on [Fe/H] was 
ever statistically significant. 

In the regime of cool giants, there is a strong evolution- 
ary effect such that metal-poor stellar populations do not 
generate M-type giants. The rich globular cluster 47 Tu- 
canae is on the cusp, containing 4 long-period variable 
stars at the tip of its giant branch at [Fe/H] k, —0.8. The 
SMC, at present-day [Fe/H] » —0.6, generates some M 
and Carbon stars, but mostly because of intermediate- 
age populations that grow very bright (and cool) asymp- 
totic giant branches. Thus, there is a sharp transition 
from excellent metallicity coverage for K giants to very 
limited metallicity leverage for M giants, exacerbated by 
the fact that M giant abundances are hard to measure. 
In M dwarfs, where stars of all metallicitics exist in our 
list, there is an interesting, strong convergence of color- 
color sequences as a function of metallicity so that G 
dwarfs have a very strong [M/H] dependence, there is 
a transition in K dwarfs, and M dwarf colors have no 
detectable [M/H] dependence. In fitting, therefore, the 
[M/H] dependence was gradually removed for cooler and 
cooler stars, for the giants because cool, metal-poor stars 
do not exist, and for the dwarfs because the [M/H] de- 
pendence removes itself empirically. 

2.4. Temperatures 

Due to our approach of fitting color-color relations in- 
ternally as a function of gravity and abundance, attach- 
ment of temperature scales could, in principle, be done 
for any color-temperature relation in any part of the pa- 
rameter space. The first iteration of this process was to 
layer color-temperature relations on top of each other 
until the whole parameter range was covered, and to 
take the median in regions where more than one rela- 
tion applied. For FG K giants lAlonso et alJ l)1999aj) and 
lAlonso et al.l l|1999bj) were used. These works include a 
specific [Fe/H] dependence, and the average of both V— I 
and V — K relations were u sed. For appropriate r uns o f 
temperatures and gravities. I VandenBerg fc Cleml (2003) 



V — I was translated to V — K via our color-co lor rela- 
tions. In a similar manner, the sy nthetic fluxes of lKuruca 
jl992b|) and lBessell et all (|1989l 119911) fluxes were c om- 
bined and translated to colors as in[Worthcv ( 1993). I n 
this case, both V — R and V — I were translated toV — K 
via the emprical color-color relations and plotted along 
with t he untweaked V — K - T e fj relations. Tokunaga 
(2000) developed average color-temperature relations for 
sequences of supergiants, gi ants, and dwa r fs usi ng lit- 
erature temperature scales. iBessell et all l)1998|) gives 
color-temperature sequences for solar-abundance dwarfs 
and giants based on different model atmospheres, and 
we also referred to the empirical co o l gian t sequences of 
iRidgwav et"al~l l)1980|) : iDvck et alJ jll996T) . The color- 
temperature sequences of iJohnsonl J19661) are also in- 
cluded. Fo r the coolest dwarfs, analysis of the data of 
iBasri et al.l 1J2000J) yielded a relation as a function of I—K 
color, specifically T eff = -460.25 x (I- K) +4323, valid 
for I -K > 2.9. 

M dwarfs with known angular diameters were also 
inclu ded in the mix. E clipsing binaries Y Y Gemino- 
rum jTorres fe Ribasl2002|) and CM Draconis IjViti et alJ 
1997) were supplemente d with inte r ferom etrically- 
derived temperatures from Bcr ger et alJ (I2006D and, in 
the ca se of Barnard's star, from iDawson fc de Robertisl 
(2004). VK ph otometry came from either our own cat- 
alog or that of iLeggett et all 112 001'). The temperature 
estimates of B erriman et alJ l)1992j) for eleven dwarfs are 
also plotted. These are more indirect temperature esti- 
mates from the ratio of the bolometric flux to the flux 
at an infrared wavelength, the total to infrared flux ratio 
method (TIRFM). The position of especially the cooler 
stars was influential in our adopting a somewhat cooler 
temperature scale around 3000 K than the bulk of the 
published calibrations. 

The fits are good to a limit ofV — K = 10.2. Since cool 
dwarfs and giants have different temperature scales, this 
corresponds to approximately T e ff = 2700 K for solar- 
metallicity giants and T e ff = 1830 K for solar-metallicity 
dwarfs. 

2.5. Bolometric Corrections 

The last item to be added was the 1^-band bolomet- 
ric correction (BC). Since they were that last item in 
the chain, BCs could be inserted as a function of color 
or of temperature and for any passband. As for tem- 
perature scales, a variety of empi rical and theoretical 
option s were intercompared. The IVandcnBc rg fc Cleml 
(2003) BCs were adopted for the middle of the temp era- 
ture range, supplemented by the lVacca et al. (1996) for- 
mula for 4.40 < log T < 4.75 for the hottest dwarfs and 
sup ergiants. At the cool end, fo r both giants and dwarfs, 
the iVandenBerg fc Cleml p003T) BCs drift slightly from 
most calibrations. For giants, we adopted the average, 
empirical-plus-theoretical iiT-band BC from Besscll et al. 
(1998), read from their Figure 20. For cool dwarfs, 
w e adopt the K - band (UKIRT IRCAM3 system) BCs 
of ILeggett et al.1 (|2001). We extended their polynomial 
slightly to reach our V — K = 10.2 cool limit. One 
subtlety regarding the Leg gett et alJ 1)200 1|) calibration 
should be mentioned. They give two polynomial fits to 
the JC-band BC, one as a function of I — K and one as 
a function of J — K. We adopt the I — K version, as 
the J — K version drifts significantly from the I — K ver- 
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Fig. 7. — Temperature-V — K calibrations for cool, solar 
abundance gian ts. L ine s are color-co ded for the calibra tions of 
lAlonso et"aU 11999a|P: IVandenBerg fc Cleml 12003^: ITokunaga I 
(2000) Tjohnsonlll96g^lWorthevlll994l) : IBessell et ali 11993) . Our 
adopted relation is shown as diamonds. 




Fig. 8. — Temperature-V — K calibrations for cool, solar 
abundance dwar fs. Lines a r e col or-coded for the calibrations of 
VandenBer g fc Cleml 120031): ITokunaga I 120001): [John son (1966); 
I Worthev I HL994): Bes sell et al.l 119981) and feasri et al.H20001) "Red 
dots with error bars are M dwarfs are from Bcrgcr et alj 120061) 
and magenta open circles are TIRF M temperatures and pho- 
tometry from |BgrrijnjMi_eL_aLl tl992|). YY Geminorum's tem- 
perature is from |Torr cs fc Ribas (2002), Barnard's star from 
Dawson & dc Robcrtis (2004), and CM Draconis's from lViti et alj 
(1997). Our adopted relation is shown as diamonds. 



sion at warmer temperatures. The cause of this drift is 
increased scatter in the J — K diagram, or, more funda- 
mentally, the fact that both J- and if-bands are on the 
red tail of the blackbody curve for the warm half of the 
temperature range covered, so that J — K as a temper- 
ature indicator has a small temperature range per unit 
error. 

Note that, plotted as a function of color, and as pre- 
dicted by synthetic fluxes, the bolometric corrections are 



Fig. 9. — V-band bolometric corrections for dwarfs and gi- 
ants. The top panel is a sequence of giants, and the bot- 
tom panel is a sequence of dwarfs, both for near-solar metal- 
licity. The symbol key is marked on the plot itsel f with 
sources Vacca et al. (1996) for hot stars, Houdashelt et al. (2000) 
for cool s t ars. IVa ndenBer g fc Cleml 120031) : IBessell et al.1 <199a) : 
IWorthev I 119941): IPlez. Brett, fc Nordlund] jl992l) : IBrettf 1199a) . 
and lLeggett et al.1 120011) , plus "PHOENI X," which refers to fluxes 
produced from the Phoenix code (Allard & Hauschildt 1995) with 
colors generated as in Worthcy ( 1994) and "Plez 1997," which 
re fers to a priva t e com munication that was subsequently published 
in Bes sell et"aT1 ^98). 



a very weak function of abundance and gravity. This is 
a degeneracy. That is, increasing a cool giant's abun- 
dance (for example) will make it redder and give it a 
larger (absolute value of the) V-band BC. Such vectors 
lie closely along the trend caused by temperature, so 
BCs are strongly covariant with T, log g, and [M/H] 
when plotted versus V — K. We exploit this for cool 
stars by adopting BCs that vary as a function of color 
alone. Gravity and abundance dependence then is in- 
herited from the gravity and abundance variations of the 
color-color diagrams. 

3. COMPARISONS AND DISCUSSION 

The wealth of comparison data that we could be check- 
ing against is too vast to illustrate completely in the 
pages of this journal, so we limit ourselves to a few key 
examples. One region of parameter space that is of keen 
interest is that of low stellar tem peratures. We check ou r 
results for coo l stars a gainst VandenBerg ^^Cleml (2003) , 
ILeieune et al.l lj!998|) . an update of the IGreer] (I1988T) 
color table used in the Yonsei-Yale isochrones l)Yi et al. I 
20011). and, for good measure, the synthetic colors of 



Worthevl IJ1994J) in Figure^] The coolest giants are very 



important for integrated-light studies of spectral features 
such as TiO that become strong only in these stars, and 
for surface brightness fluctuation (SBF) magnitudes, es- 
pecially at red colors, that depend on these stars because 
of the L 2 dependence of an SBF magnitude. 

In Figure ITUI our fits are shown as black lines. They 
were fitted to the B - V - V - K and V - I - 
V — K diagrams, so it is no surprise that they still 
fit in this color-color plane. The updated-Green cali- 
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Fig. 10. — The B — V, V — I color-color diagram for M stars. 
Red dots are stars with [Fc/H] > —0.2 and blue dots are stars 
with [Fc/H] < — 1.2 with intermediate stars cyan. Giants are open 
symbols, which dwarfs are filled. The data are unculled. Calibra- 
tions for giant and dwarf color-color sequences are drawn in solid 
for [Fe/H] = 0, dashed for [Fe/H] = -1, and dotted for [Fe/H] 
= —2. The c olor codes for different authors are noted in the figure 
("Leje une"is Lcicunc ct al. 11998 ), "Green" is the upda ted IGreenl 
tl98S) table, and "V&C" is h/andenBerg fc Cleml J200^ 1. 



Fig. 11. — R — I is plotted as a function of V — K. Stars of dif- 
ferent giant/dwarf status and abundance are plotted with different 
symbols according to the key. Stars with [Fe/H] > —0.2 are con- 
sidered metal rich, stars with [Fe/H] < —1.2 are considered metal 
poor, and stars between these values are c ons idered intermediate 
in metallicity. Lines are coded as in Figure Hoi This color was not 
fit during the calibration process. 



bration follows an extrapolation of the giant sequence 
off into regions not occupied by stars, while the dwarf 
sequence for solar abundance follows the stars very 
well. There is considerable metallicity dependence in 
the Green calibration that the stars do not appear to 
share. The VandenBerc & Clem ( 2003) sequences follow 
both dwarfs and giants fairly well, with a fairly good 
(small) abundance dependence. The oscillations in the 
solar metallicity giant track are a reflection of actual 
values in thei r data tab les. The coolest temperature 
reached by iVandenBerg fc Cleml l|2003[) is 3000 K. The 
iLeieune et al. ( 1998) calibration is based on corrected 
synthetic fluxes. In this case, the dwarfs and giants track 
together with little or no gravity separation until, at a 
temperature well within the tabulated range of applica- 
bility, the values become wild. 

Besides author comparisons, another way to check our 
results is to plot colors that were not fitted explicitly to 
see if the implicit color dependences are correctly mod- 
eled. R — I is one such, and is illustrated in Figure ^] 
For this color, the fits were versus V — R and V — I, 
for slightly different samples of stars. The R — I fitted 
tracks fall among the stars fairly well, except for a rever- 
sal seen around V — K = 1.5 where the giants become 
w 0.02 mag redder than the dwarfs. This 0.02 mag shift 
is probably incorrect, but it gives a valuable indication 
of the reliability of the color-color fits. 

Reliability must be a function of temperature regime. 
One particular troublesome area is that of K-type dwarfs 
and the damping of the abundance sensitivity going to- 
ward cool stars. In Figure IT2l a T e ff of 5000 K corre- 
sponds to V — K rs 2.2 and T e // = 4000 K corresponds 
to V — K = 3.4. It is clear that the magic combination of 
full photometry plus a good abundance estimate is lack- 



ing from our data set for dwarfs. Note the general lack of 
open symbols (dwarfs) compared to filled (giants). This 
lack of data means that the metallicity dependence of 
the dwarfs is inherited from the plentiful giants in this 
temperature regime; an undesirable feature. At redder 
colors, as their surface gravities diverge, the dwarfs and 
giants separate in color. Simultaneously, the metallicity 
dependence appears to reverse, at least in the giants. It 
is unclear from the present data what should be happen- 
ing with the dwarfs. The polynomials do their best to 
smoothly flow through all of this, but we judge it un- 
likely that they have truly captured the essence of the 
color behavior in this regime, as it is not clear to our 
eyes exactly what should be happening (it seems likely 
that some of the photometry is bad). This particular 
problem is an issue only for U — B and B — V colors, 
although the paucity of K dwarf data is of some concern 
for all colors, as the metallicity dependence is relatively 
unconstrained. 

The principle source of error in the color-color fits is 
finding a suitable polynomial to follow the various twists 
and turns that the colors take. We fit the colors in five 
segments, with multiply-redundant overlap in color, and 
used the overlap regions to estimate the error from poly- 
nomial fitting. With typically hundreds of stars available 
for each fit, random photometric uncertainty is not a con- 
cern (though, of course, systematic uncertainty is). The 
median fit uncertainty over all temperatures, gravities, 
and abundances is listed in Tabled We also thought it 
useful to propagate errors in the final subroutine so that 
uncertainties in the effective temperature scale could be 
translated to uncertainties in color. For this we used the 
various T e ff relations plotted in Figures and |H1 and a 
couple of others to roughly estimate a percentage error 
as a function of temperature. This is given in Table [21 
Note that the errors in Table|2]for cool stars are more ap- 
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TABLE 2 
Temperature 
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Assumed 



Fig. 12. — The figure show a small section of the B — V, 
V — K color-color diagram. Dwarfs are drawn with larger sym- 
bols than giants for emphasis according to the key in the figure. 
Stars with [Fe/H] > —0.2 are considered metal rich, stars with 
[Fe/H] < —1.2 are considered metal poor, and stars between these 
values are considered intermediate in metallicity. Calibrations for 
giant and dwarf color-color sequences are drawn in solid for [Fc/H] 
= 0, dashed for [Fe/H] = -1, and dotted for [Fe/H] = -2. At red 
color, the dwarfs follow the bluer B — V tracks. This is a region of 
uncertainty, as discussed in the text. 



TABLE 1 

Median 
Polynomial Fit 

Uncertainty 



Color 


a (mag) 


U-B 


0.071 


B-V 


0.017 


V -R 


0.010 


V -I 


0.011 


J -K 


0.004 


H - K 


0.002 



plicable to giants than dwarfs; dwarf temperatures seem 
more uncertain than those of giants, but we didn't have 
enough dwarf calibrations to estimate this very well, so 
we left it alone. For color / with color error 07 and a 
temperature error or, errors propagate in the elemen- 
tary way: 






(1) 



Color-color diagrams have been used to derive a "color 
excess" fro m which can be infe rred a value for the dust 
extinction (Morga n et al.lll953j) . The metallicity depen- 
dent color-color fits of this paper offer a general, if not 
overly precise, method of generating a color-color plot for 
any color combination as a function of abundance and 
gravity. The classic U — B, B — V diagram is shown in 
figure IT^l for dwarfs only. The double inflection redward 
of zero color represents the rise and fall of the Balmer 
break in B-type through A- and F-type stars. A defect 
of this method is that it only works on clusters that have 
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A-type stars, that is, ones younger than about 1 Gyr that 
still have dwarfs that hot. Interestingly, there is an addi- 
tional color inflection in the M dwarfs, roughly between 
4000 K and 3000 K, that may allow independent redden- 
ing estimates for old clusters that have deep photometry. 
This inflection exists in almost every color, although the 
U band presents the most dramatic manifestation of it. 
If the U band is utilized, the coolest stars involved 
have U — I = 5.65 mag accordi ng to our colo rs and 
Mi = 9.0 according to diagrams in lLegget t (1992). This 
leads to an absolute Mu — 14.65 mag. If a modest tele- 
sc ope can reach U = 23 as the KPNO 2.1-meter did 
in IKaluznv fc Rucinskil IJ1995J) . then the U flux is read- 
ily detectable to about 500 pc distance. For reference, 
the nearest ancient open cluster is M67 at about 800 pc, 
so one would need a slightly larger telescope or better U 
sensitivity for U to be useful. However, redder colors can 
also be made to work at about the same confidence level 
relative to the fitting errors. The fitting errors are shown 
as error bars on points in Fig. 1131 as is another color se- 
quence shifted by Ay = 0.1, shown as a line. We judge 
that this extinction is the smallest that can be detected 
at all simply using the color-color fits we present, and so 
is not particularly competitive with other methods as it 
stands. Interestingly, at redder passbands, the M-type 
deflection becomes less pronounced but the errors also 
decrease so that the Ay = 0.1 extinction vector stays 
at about the same statistical significance in most color- 
color planes. This does not solve the problem, however, 
because observational measurement error becomes larger 
than the fit error at JHK wavelengths. Future refine- 
ments to this reddening estimation method are possible 
and should be encouraged. 

4. SUMMARY AND CONCLUSION 

Johnson/Cousins photometry was combined with lit- 
erature [Fe/H] estimates to fit color-color diagrams as a 
function of gravity and abundance. Literature-average 
temperature and bolometric correction scales are at- 
tached to provide a global color-temperature relation for 
stars with -1.06 < V - K < 10.2. The RI magnitudes 
are in the Cousins system, and JHK magnitudes are in 
the Bessell homogenized system. The complete color- 
temperature table and a Fortan interpolation program is 
available at http://astro.wsu.edu/models/ 

Several areas of improvement were noted in the main 
body of the paper, including filling photometry gaps, ob- 
taining more accurate and on-system photometry, know- 
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ing better log g and [Fe/H] values, improving the statis- 
tics for data-impoverished groups of stars such as K 
dwarfs, applying small tweaks in the processing pipeline, 
and obtaining better empirical temperature and bolo- 
metric correction relations, especially for supergiants and 
M stars. 

A way to estimate dust extinction from M dwarf colors 
arises from an inflection that exists in most colors rela- 
tive to V — K. Unlike the classic UBV method, it can be 
used in old star clusters, but it does not seem to promise 
much, if any, increase in accuracy for clusters where both 
methods apply. The most sensitive band relative to pho- 
tometric error for the new extinction measure is the U 
band, but if the U band is employed then clusters must 
be within a few hundred parsecs for ground-based obser- 
vatories to able to measure adequate U fluxes. 



Fig. 13. — The final color-color calibration for dwarfs only is 
shown in the U — B, B — V plane as dots with error bars attached, 
where the error bars include fitting uncertainties and T e jj uncer- 
tainties. An additional line is shown that represents the color shift 
due to dust screening of Ay =0.1 mag. For illustrative purposes, 
a vector for Ay = 1.0 mag is also sketched. The approximate 
blue limits of the bluest dwarfs at the main sequence turnoffs of 
isochrones of various ages are marked with the corresponding ages. 
Near the red bump feature, stellar effective temperatures (degrees 
K) are indicated. 
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